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ABSTRACT 

Observational evidence is presented for periodically variable irradiation of secondary compo- 
nents. This results in strongly modulated mass outflow. Superhumps are then due to enhanced dissi- 
pation of the kinetic energy of the stream. 

Qualitative interpretation of superhump periods and their variations is also presented. 
Key words: accretion, accretion disJcs - binaries: cataclysmic variables, stars: dwarf novae 

1. Introduction 

Superhumps are common among dwarf novae of the SU UMa subtype during 
their superoubursts (Warner 1995, Hellier 2001). They are also observed in the case 
of the so-called permanent superhumpers (cf. Patterson 1999). The superhump pe- 
riods are slightly longer than the orbital periods and their amplitudes are - typically 
- 0.3 mag. 

The commonly accepted "tidal" model explains superhumps as being due to 
tidal effects in the outer parts of accretion disks leading - via the 3:1 resonance - 
to the formation of an eccentric outer ring undergoing apsidal motion. This model 
and, in particular, the results of numerous 2D and 3D SPH simulations (cf. Smith 
et al. 2007 and references therein) successfully reproduce the observed superhump 
periods and correlations of the superhump period excess with the orbital period and 
the mass-ratio. Moreover, those correlations find natural interpretation in the con- 
text of structure and evolution of the secondary components (Patterson 1998,2001, 
Patterson et al. 2005). 

On the other hand, however, the tidal model fails to reproduce the amplitudes 
of superhumps: The numerical SPH simulations produce "superhumps" which 
have amplitudes ~ 10 times too low, compared to the observed amplitudes (Smak 
2009a). This implies that the nature of superhumps requires another explanation. 



2 



A. A. 



Recent analysis of superhumps in Z Cha (Smak 2009b) showed that - contrary 
to the "tidal" interpretation - they are due to strongly modulated mass transfer rate. 
The question then appears: what is the origin of those modulations? The answer to 
this question is given in the present paper. 

In Section 2 we present observational evidence for variable irradiation of sec- 
ondary components. Its consequences are discussed in Section 3 in terms of the 
irradiation modulated mass outflow. In Section 4 we discuss the relationship be- 
tween the irradiation periods and the observed superhump periods and present gen- 
eral discussion of their variability. Results are summarized in Section 5. 



Many years ago Schoembs (1986) discovered that the luminosity of OY Car 

during its November 1980 superoutburst varied with period close to the beat pe- 
riod. Until now, however, the significance of this discovery has not been fully 
appreciated. 

The see whether such a behavior is common among other dwarf novae a search 
was made through the literature for light curves well covering all beat phases. Fif- 
teen such cases were found, including 12 dwarf novae observed during their super- 
outbursts and 3 permanent positive superhumpers. They are listed in Table 1, where 
the first six stars, which are eclipsing systems, are arranged in order of decreasing 
orbital inclinations, while the remaining, non-eclipsing systems - in alphabetical 
order. For eight of them it was possible to use the pure disk magnitudes measured 
away from eclipses and from superhumps. For the others - only nightly mean mag- 
nitudes were available. 

We analyze those light curves in the following way. To begin with, using the 
orbital and superhump elements, we calculate the beat phases 



In three cases, where no orbital periods were available, their values had to be de- 
termined from their superhump periods using formula given by Menninckent et al. 
(1999). Consequently an arbitrary zero-point of i^t> had to be adopted. 

Anticipating further considerations it may be worth recalling that the beat phase 
refers to the beat period which is defined by the other two periods as 



2. The Missing Link 



(1) 



1 



1 



1 



Pb 



Porb 



Psh 



(2) 



The observed light curves are fitted with a simple formula: 



m = m{to) + — {t-to) - A cos {<Sfb-<Sfb^) , 



(3) 
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Table 1 
Photometric Data 



Star 


i 


Data source 




lYUMa 


86.8 


Patterson et al. (2000a), Fig.2, JD245 1562-569 


2 


DVUMa 


84.0 


Patterson et al. (2000b), Table 1 + Fig.3, 1997 


3 


OYCar 


83.3 


Schoembs (1986), Fig.l 


2 


ZCha 


80.2 


Kuulkers et al. (1991), Table 3 


3 


WZSge 


75.9 


Patterson et al. (2002), Table 1, JD2452126-137 


3 


UGem 


69.0 


Smak and Waagen (2004), JD 2446344-376 


3,4 


V603 Aql 




Patterson et al. (1993a), Table 1 


3 


VY Aqr 




Patterson et al. (1993b), Fig.7, JD2446558-566 


2,5 


TTAri 




Wuetal. (2002), Fig.l 


2 


TTBoo 




Olech et al. (2004), Fig.3 


2,5 


V503 Cyg 




Harvey et al. (1995), Table 3, JD2449601-609 


3 


VWHyi 




Schoembs and Vogt (1980), Fig.2a, JD2443808-815 


2 


BKLyn 




Skillman and Patterson (1993), Table 1, 1992/93 


1,3 


TU Mon 




Stolz and Schoembs (1984), Fig.2 


2 


KSUMa 




Olech et al. (2003), Fig.3 


2,5 



Remarks: 1. PG 0917+342, 2. pure disk magnitudes, 3. mean magnitudes, 4. mean 
points in 0.1 phase bins, 5. zero-point of arbitrary. 



where the first two terms describe the slow brightness decline during superoutburst, 
while the third - the expected modulation with the beat period. The residual mag- 
nitudes are then determined as 



Am 



m 



(4) 



Results are Usted in Table 2 and shown in Figs.l and 2. As we can see there 
is an obvious dependence on the orbital inclination: (1) There are only three cases 
showing large variations of Am with (|)^, with full amplitudes 2A ~ 0.35 — 0.40 
mag. Those are the deep eclipsers with highest orbital incUnations (OY Car, lY 
UMa and DV UMa). Worth noting is that they have also very similar phases of 
maximum. (2) The other three eclipsing systems with lower inclinations (Z Cha, 
WZ Sge and U Gem), show much smaller amplitudes which, in view of their large 
errors and large scatter of points seen in Fig.l, can hardly be considered significant. 
Plotted in Fig.2 is also another eclipsing dwarf nova XZ Eri, which was analyzed 
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earlier by Uemura et al. (2004). During one of its superoutbursts they detected low 
amplitude variations with a period P = 4.7 d, but no measurable variations with the 
beat period Pb = 2.33d. It is plotted in Fig.2 with A = at ? = 80.2 as an asterisk. 
(3) The non-eclipsing systems, with low incUnations, do not show any significant 
variations. 



Table 2 
Results 



Star 


A (mag) 


xmax 


Star 


A (mag) 


^.max 


lY UMa 


0.1 8 ±0.04 


0.67 ±0.04 


VY Aqr 


0.02 ±0.05 




DVUMa 


0.17±0.06 


0.60 ±0.06 


TTAri 


0.02 ±0.02 


0.38±0.12 


OYCar 


0.20±0.07 


0.65 ±0.07 


TTBoo 


0.05 ±0.03 




ZCha 


0.06 ±0.09 


0.28 ±0.25 


V503 Cyg 


0.02 ±0.06 


0.10±0.21 


WZSge 


0.03 ±0.03 


0.49 ±0.1 3 


VWHyi 


0.05 ±0.03 


0.20 ±0.20 


UGem 


0.04 ±0.04 


0.35±0.18 


BKLyn 


0.02±0.01 


0.84±0.17 








TUMon 


0.01 ±0.08 


0.59 ±0.23 


V603 Aql 


0.02 ±0.02 


0.35±0.19 


KSUMa 


0.02±0.01 





There appears to be a simple interpretation of those results. Disks in dwarf 
novae during their superoutbursts and in permanent superhumpers have strongly 
non-axisymmetric structure, rotating - in the observer's frame - with a period Prot , 
which is related, but only approximately equal (see Section 4), to the observed beat 
period Pb- In particular, the effective geometrical thickness of the disk z/r must 
depend on the position angle or - in the observer's non-rotating frame - on the 
phase of Prot ■ If so, the observed periodic light variations with period Pb can be 
interpreted as being due to variable obscuration of the central parts of the disk, 
resulting from variable z/r. 

It is obvious that such variations can be seen only at orbital inclinations fulfill- 
ing the condition: cos / « z/r. As can be seen from Fig.2, the three systems show- 
ing such variations have orbital incUnations corresponding to z/r = 0.06 — 0.12. 
This is, indeed, typical for hot disks with M = 10^^ - lO^^g/s (cf. Smak 1992). In 
the case of Z Cha and XZ Eri, however, we get z/r = 0.17 which would require 
unacceptably high accretion rate M = lO^^g/s. 

Modulations of disk luminosities with P^ot (in the non-rotating frame of the 
observer) imply variable irradiation of the secondary component (in its rotating 
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Fig. 1. Residual magnitudes versus beat phase for stars listed in Table 2. Solid lines are cosine 
curves with A and obtained from solutions with Eq.(3). 
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60 70 80 

inclination 



Fig. 2. The amplitude as a function of inclination for stars listed in Table 2. Non-eclipsing systems 
with unknown inclinations are plotted between i = 50° and 60° . An asterisk denotes XZ Eri. Marked 
above are values of z/r = cos (' . 

frame) with "irradiation period" P,>r wiiicii - in analogy to Eq.(2) - is defined by 

1 _ 1 1 

Prot ^orb Pirr 

and therefore is approximately (see Section 4) equal to the observed superhump 
period. 



3. Irradiation Modulated Mass Outflow 



Irradiation of the secondary component was discussed by many authors as an 
important factor controlling the mass outflow from the inner Lagrangian point L i . 
The main problem here is connected with the fact that the equatorial parts of the 
secondary, including Li , are in the shadow cast by the disk and therefore are not 
directly affected by irradiation. Results of very crude model calculations (Smak 
2004) showed that (1) the hot material flowing from irradiated regions can reach 
the vicinity of L i , and (2) its temperature at L i is still high enough to produce 
substantial enhancement in the mass outflow rate. More recent model calculations 
by Viallet and Hameury (2007) confirmed the first conclusion, but not the second: 
according to them the flow, when it reaches L i , is already too cool to produce any 
significant effects. It appears, however, that this result was mainly due to one of 
their assumptions. The rate of cooling of an isothermal layer is determined by its 
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surface effective temperature Te which - due to that coohng - is lower than the 
temperature T of the layer itself. Viallet and Hameury assumed Te = T thereby 
overestiming the cooling effects. Further model calculations are needed. 

The superhump phase of maximum irradiation of the secondary component 
^nmx related to the observed phase of maximum luminosity (sf^ (see previous 
Section) by (|)"|''^ = — (|)™^ . This means that the maximum irradiation preceeds the 
next superhump by 

A/ = <SfTPsh = ^r^Porb. (6) 

"orb 

This delay consists of two terms: 



At = At flow + Atstr , (7) 

where Afj/ow is the time needed for the flow to reach Li, and At str is the time 
needed for the stream to reach the point of impact. 

The last equation can be used to estimate At flow From lY UMa, DV UMa and 
OY Car, using their average values of < (j)^"^ >= 0.64 and < Psh/Porb >= 1-028, 
we get At = 0.66 Porb- Values of Atgfr can be obtained from stream trajectory 
calculations. With < ^ >= 0.13 we get from At str = 0.19 Porb at the edge of 
the disk to Atstr = 0.23 Pgrb at a point half way between disk edge and its center. 
Inserting those values into Eq.(7) we obtain: AtfUm = (0.43 — 0.47) Porb- 

Two comments are worth to make at this point. First, that surprisingly similar 
values: Atfi ow ~ 0.4 Porb were predicted by our crude model calculations (Smak 
2004, Table 2). Secondly, that - should the flow time be actually much longer 
- then our result would have to be rewritten in a more general form as Atfiow = 

{0.43-0.41} Porb + nPorb- 



4. The Superhump Periods and their Variations 

In the previous Sections we discussed only the effects of periodic variations of 
z/r which occur on a short time scale corresponding to Psh- Let us now consider 
its variations on a longer time scale. During superoutbursts the mass transfer and 
accretion rates change considerably resulting in significant variations of the irradi- 
ating flux and the effective geometrical thickness of the disk (averaged over all po- 
sition angles). Generally we have: Firr ~ M and z/r M (cf. Smak 1992). There- 
fore: dFirr/dM > and d{z/r)/dM > 0. The flow time depends on both those 
factors and, in particular, we expect dAt flow /^firr < and dAtfiow/^{z/r) > 0. 
The flow time is then predicted to vary as 



dAt 



flow 



dt 



dAtfiow d(z/r) ^ dAt flow dFi 
d{z/r) dM 



+ 



dFirr dM 



dM 
dt 



(8) 
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We note that the first term in the square brackets (with z/r) is positive while the 
second term (with Firr) is negative. Therefore, depending on which of the two 
factors is dominant, the expression in the square brackets can be either positive or 
negative. 

Taking into account variations of Atfja^-, and assuming that At^tr is practically 
constant, we obtain the following equation connecting the observed superhump 
period with irradiation period (introduced in Section 2 via Eq.5) 



fsh — 



dAtfiow 

1 7 

dt 



(9) 



As we can see, the two periods are equal only when At flow = const (in such a case 
we also have Pb = Pm)- When ^!Ufiow changes with time the two periods differ by 
a small amount, depending on the value and sign of dAtfiow/dt. 

Let us now turn to the problem of superhump period variations. Up to mid 
1990-ies nearly all dwarf novae observed until that time showed superhumps with 
periods which were decreasing during superoutburst (see Warner 1995). This was 
interpreted (cf. Lubow 1992) as being due to decreasing disk radius. Such an ef- 
fect was indeed predicted by the thermal-tidal instability model of Osaki (1989). 
In particular, his calculations showed that the contraction of the disk begins at the 
very beginning of the superoutburst and continues throughout its main part. This 
interpretation, however, must be abandoned for two reasons. First, there is no obser- 
vational evidence for such disk radius variations. Secondly, superoutbursts are not 
due to thermal-tidal instability, but due to a major enhancement in the mass transfer 
rate (Smak 2008). As a result, the accretion is quasi-stationary, with strongly en- 
hanced accretion rate, causing disk radius to be close to the tidal radius and remain 
constant throughout the entire superoutburst. 

Situation became more complicated with the discovery of a group of dwarf 
novae with shortest orbital periods showing increasing superhump periods (Kato et 
al. 2001,2003, Olech et al. 2003, and references therein) and even more so with the 
discovery of several examples of more complex alternating period variations (e.g. 
Olech et al. 2004, Rutkowski et al. 2007). 

The new interpretation of superhumps emerging from considerations presented 
above adds another dimension to this problem. From Eq.(9) we get 

dPsh _ Psh dPjrr d^Atflow 
dt ~ Pirr dt ^ dt^ ■ ^ 

This equation shows that superhump periods can change due to (1) variations of 
Pirr (resulting from variations of Prot), and (2) non-linear variations of i!Ufiow (as 
described by its second derivative). It explains why the observed period variations 
can be so complex. 

At the present moment we can present - mainly for illustrative purposes - only 
some crude, quaUtative considerations. For this purpose let us arbitrarily assume 
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(1) that Pirr is constant, and (2) that the expression in the square brackets in Eq.(8) 
remains constant or changes with time only slightly so that the main contribution 
to the second derivative of A///ow in Eq.(lO) comes from variations of M. If so, 
Eq.(lO) can be replaced with 



dt ~ 



dAtfiow d{z/r) ^ dAtfiow dFir 



d{z/r) dM dFirr dM 



d^M 



The accretion rate decreases during the main part of the superoutburst. In particular, 
we have dlogM/dt Ri const (cf. Fig. 3 in Smak 2008) which implies d^M/dt^ > 
0. The resulting variations of Psh, however, will be either positive or negative, 
depending on the sign of the expression in the square brackets: When z/r is the 
dominant factor we predict dPsh/dt > 0. Conversely, when Fjrr is the dominant 
factor we expect dPgh/dt < 0. 

Turning to observations we recall that superhump periods of dwarf novae with 
longer periods are generally decreasing, while those of dwarf novae with shortest 
periods are increasing (Kato et al. 2003, Fig. 14, Olech et al. 2003, Fig.9). The 
transition occurs around Ptr ^ 0.06 d. Our predictions could then suggest that the 
dominant factor in the first case is z/r , while in the second - Fjrr . In addition, there 
is a group of dwarf novae showing alternating variations of Pgh (Olech et al. 2003, 
2004, Rutkowski et al. 2007). Their periods are close to Ptr and this could suggest 
that in their case z/r nd Firr are equally important, resulting in a more complex 
behavior. 

Needless to say, detailed model calculations will be needed to confirm this 
simple interpretation and to perform the quantitative analysis of the problem. 



5. Discussion: the Nature of Superhumps 



Combining earlier evidence showing that superhumps in Z Cha are due to 
strongly modulated mass transfer rate (Smak 2009b) with evidence presented in 
the present paper we propose the following new intepretation of superhumps: 

(1) Irradiation of secondary components is observed to vary periodically with 
period Pirr which is close (but not exactly equal) to the observed superhump period 
Psh- The basic "clock" mechanism must then be the same as in the case of the 
"tidal" model. 

(2) Variable irradiation of the secondary component results in mass onflow be- 
ing strongly modulated with period equal to the observed superhump period. 

(3) Superhumps are then due to enhanced dissipation of the kinetic energy of 
the stream. 

Two kinds of models are now needed for a more meaningful, quantitative inter- 
pretation of superhumps and, in particular, of their variable periods: 
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(1) Realistic 3D disk models, including the effects of stream overflow and vari- 
able mass transfer rate. They are needed (?) to identify the nature of the "clock" and 
(ii) to determine z/r and Firr as functions of time. 

(2) Models describing the irradiation modulated mass outflow. They are needed 
(0 to determine Atfiow and M as functions of z/r and Fjrr (see Section 4) and (//) 
to study the relationship between P,>r and Pgh ■ 

It can be hoped that once such models become available it will be possible 
to discuss not only the common superhumps and their periods, but also the other 
members of the "superhump zoo" (cf. Patterson et al. 2002). 

Acknowledgements. The author is grateful to Dr. Arkadiusz Olech for helpful 
discussions. 
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